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ABSTRACT 

We present high-precision radial velocity (RV) observations of four solar-type (F7- 
G5) stars - HD 79498, HD 155358, HD 197037, and HD 220773 - taken as part of 
the McDonald Observatory Planet Search Program. For each of these stars, we see 
evidence of Keplerian motion caused by the presence of one or more gas giant plan- 
ets in long-period orbits. We derive orbital parameters for each system, and note the 
properties (composition, activity, etc.) of the host stars. While we have previously 
announced the two-gas-giant HD 155358 system, we now report a shorter period for 
planet c. This new period is consistent with the planets being trapped in mutual 2:1 
mean-motion resonance. We therefore perform an in-depth stability analysis, placing 
additional constraints on the orbital parameters of the planets. These results demon- 
strate the excellent long-term RV stability of the spectrometers on both the Harlan J. 
Smith 2.7 m telescope and the Hobby-Eberly telescope. 

Subject headings: planetary systems — planets and satellites: dynamical evolution 
and stability — stars: individual (HD 79498, HD 155358, HD 197037, HD 220773) — 
techniques: radial velocities — astrobiology 



1. Introduction 



Beginning with the first radial velocity exoplanet detections ( Latham et al.lll989l : lMavor Sz Queloz 



1995), exoplanet surveys have identified a large number of gas g iant planets on very short orbits 



(see lUdrv &: SantosI 120071 : iBorucki et al.ll201ll : IWright et al.ll201ll . for complete details). However, 
to answer the fundamental question "how common are planetary systems analogous to our solar 
system?" requires years of data from dedicated exoplanet surveys to achieve the observational time 
baseline necessary to detect the true Jupiter analogs-those giant planets which do not undergo 
significant inward migration after the dissipation of the protoplanetary disk. Planets with periods 
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greater than approximately 2 years are particularly valuable, as they ar e outside the detection lim- 



its of space-based transit searches such as Kepler (iBorucki et alj 120101 ) and CoRoT (|Baglin et al 
20031 ). whose mission timescales are too short to confirm such detections. 



Identifying long-period Jovian planets is essential to constraining theor etical and observationa l 
results on the true distribution of giant planets versus orbital separation. ICumming et al.l (j2008l ) 
find appro xmiately 10.5% of FGK s tars host gas giant planets with periods between 2 and 2000 
days, and IWittenmver et al.l (|2011al ) claim a 3.3% occurrence rate for "Jupiter analogs," which 
they define as gia nt planets with e < 0.2 and P > 8 years. While microlensing surveys (e.g. 
Gould et al.ll2010l ) derive a higher giant planet fraction, their results are roughly consistent with 
RV surveys if we consider the differences in sensitivity between the two methods. Microlensing 
planet searches are particula rly sensitive to p lanets near the Einstein radius (ty pically 2-4 AU) 



wher e core accretion models (iBenz et al.ll2008l ) and extrapolations of RV results (ICumming et al 



20081 ) agree there should exist a significant population of Jovian planets. The expansion of the time 
baselines of RV surveys to include these long-period planets is essential to reconcile the statistics 
of RV and microlensing planet searches. 

A more complete census of long-period gas giants also places strong constraints on theories 
of planet formation and migration. The frequen cy of Jovian planets offers information on the 
efficiency of core accretion in protoplanetary disks (IDodson-Robinson et al.ll2009l : iMann et al.ll2010l : 
Bromley Kenyonll201ll ). and the mass-period distribution will test pred ictions that planets with 
m sin i > Mj do not mig rate as far as Neptune- or Saturn-mass planets (jPel Popoulo et al.ll2005l : 
Bromley KenyonI l201ll ) . Additionally, enlarging the sample size of giant planets allows for a 
more robust examination of correlations between the properties of exoplanetary s ystems (e.g. planet 
frequ ency, m ass, eccentricity) and those of their stellar host s, such as metallicity (jFischer &: Valenti 
20051 ). mass (j Johnson et al.ll201ll ). and galactic dynamics (jEcuvillon et al.l 120071 ). 



The McDonald Observatory Planetary Search (jCochran &: Hatzeslll993l ) has been conducting 
a high-precision radial velocity survey to identify substellar companions around FGK stars with 
the Harlan J. Smith 2.7 m telescope since 1987. Since our migration to a cross-dispersed echelle 
spectrograph in combination with an I2 absorption cell in 1998, and with the addition of the Hobby- 
Eberly Telescope in 2001, we have successfully monitored hundreds of stars with velocity precision 
of ~ 3 m/s, giving us a 14 year observational time baseline. As a result, we now see evidence of 
planets with periods of 5 years or more, a demographic that is still underrepresented in exoplanet 
discoveries, accounting only 38 of the 518 planets listed in the exoplanets.org database (as of 5 
January 2012). In this paper, we present two such objects-HD 79498b and HD 220773b-a shorter- 
period Jovian planet around HD 197037, and updated orbital parameters for the two-planet system 
surrounding HD 155358. 



2. Sample and Observations 

The McDonald Observatory Planet Search Program currently monitors over 250 (mostly FGK) 
stars with the 2.7 m Harlan J. Smith Telescope for RV variations due to planetary companions. 
The survey is magnitude limited to V ^ 10, and regularly achieves RV precision of 3-6 m/s. The 
large number of nights provided for our program results in excellent temporal coverage and long 
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observational time baselines for our targets; the objects presented here each have 30-120 data points 
over 7-10 years. Thus, our sensitivity is more than sufficient to detect RV signatures comparable 
to that of Jupiter, which would require observations at our level of precision over 12 years. 



2.1. 2.7 m Telescope Observations 



The velocities for HP 79498 and HD 197037 were obtained with the Smith Telescope's Tull 
Coude Spectrograph (jTull et al.lll995l ) using a 1.8" slit, giving a resolving power of R = 60,000. 
Before starlight enters the spectrograph, it passes through an absorption cell containing iodine 
(I2) vapor maintained at 50° C, resulting in a dense forest of molecular absorption lines over our 
stellar spectra from 5000-6400 A. These absorption lines serve as a wavelength metric, allowing 
us to simultaneously model the instrumental profile (IP) and radial velocity at the time of the 
observation. For each star, we have at least one high-S/N iodine-free template spectrum, which we 
have deconvolved from the IP using the Maximum Entropy Method, and against which the shifts 
due to the star's velocity and the time-variant IP are modeled. Our reported RVs are measured 
relative to this stellar template, and are further corrected to remove the velocity of the observatory 
around the solar system bar ycenter. All radial velocities have been extracted with our pipeline 
AUSTRAL (jEndl et al.ll200d ). which handles the modeling of both the IP and stellar velocity shift. 



2.2. HET Observations 



O ur RV data for H P 155358 and HD 220773 were taken with the High Resolutio n Spectrometer 



(HRS : ITu11 et al.lll998l ) on the queue-scheduled 9.2m Hobby-Eberly Telescope (HET: lRamsev et al 
19981 ). As with the 2.7 m observations, the HET/HRS spectra are taken at R = 60,000 with an I2 



absorption cell. The fiber-fed HRS is located below the telescope in a temperature-controlled room. 
Separate stellar templ ates were obtained wi th HRS for these stars. Details of our HET observing 
procedure are given in lCochran et all ( 2004 ). 



While we use HRS spectra exclusively for obtaining RVs for HD 155358 and HD 220773, we 
have l2-free spectra from the 2.7 m telescope for these stars as well. These spectra serve two 
purposes. First, they allow us to determine the stellar parameters for all five stars using the same 
instrumental setup. Also, the Tull 2.7 m coude spectrograph provides Ca H and K indices, which 
contain information as to the activity levels of the stars. 

Tables [TJH] list the relative velocities and their associated uncertainties for HD 79498, HD 
155358, HD 1 9 7037 and HD 220773, respectively. Table [2] includes observations published in 
Cochran et al.l (j2007l ). but since all of our spectra have been re-reduced with our most up-to-date 
methods, the velocities presented here have a higher precision. 
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3. Analysis and Orbit Modeling 

3.1. Host Star Characterization 

We determine the stellar parameters of our targets according to the procedure desc ribed in 
Brugamver et alj (|201lh . The method relies on a grid of ATLAS9 model atmospheres (jKurucz 
19931 ) in combinati on with the local thermodynamic (LTE) line analysis and spectral synthesis 
program MOOgEI (Sneden 1973 ). Using the measured equivalent widths of 53 neutral iron lines 
and 13 singly-ionized iron lines, MOOG force-fits elemental abundances to match the measured 
equivalent widths according to built-in atomic line behavior. Stellar effective temperature is deter- 
mined by removing any trends in equivalent widths versus excitation potential, assuming excitation 
equilibrium. Similarly, we compute the stellar microturbulent velocity by eliminating trends with 
reduced equivalent width (= Wa/A). Finally, by assuming ionization equilibrium, we constrain 
stellar surface gravity by forcing the abundances derived from Fe I and Fe II lines to match. 

We begin our stellar analysis by measuring a solar spectrum taken during daylight (through 
a solar port) with the same instrumental configuration used to observe our targets. The above 
procedure yields values of T^// = 5755 ± 70 K, log^ = 4.48 ± 0.09 dex, ^ = 1.07 ± 0.06 km/s, 
and loge(Fe) = 7.53 zt 0.05 dex. We then repeat this analysis for our target stars, using the l2-free 
template spectra. We note that our derived metallicities are, as is conventional, differential to solar. 
The high S/N and spectral resolution of our stellar templates allow us to make robust estimates of 
each star's effective temperature, log 5, metallicity, and microturbulent velocity, which we include in 
Table [5l We also include photometry, p arallax data, and spectral types from the ASCC-2.5 catalog 



(Versi on 3. iKharchenko &: Roeseijl2009l ). as well as age and mass estimates from lCasagrande et al 
(2011 



In addition to basic stellar parameters, we use the standard Ca H and K indices to evaluate 
the hypothesis that our observed RV signals are actually due to stellar activity. For HD 79498 and 
HD 197037 we have time-series measurements of the Mount Wilson Shk index. For each of these 
stars, we obtain Shk simultaneously with radial velocity, so we inc lude those values in the RV 
tables. From the average value of these measurements, we derive (via lNoves et al.lll984l ) logi?^^, 
the ratio of Ca H and K emission to the integrated luminosity of the star, which we include in 
Table m The activity indicators are discussed in detail for each star below, but we note here that 
the sample overall appears to be very quiet, and we have no reason to suspect stellar variability as 
the cause of the observed signals. 

As a second fail-safe against photospheric activity mimicking Keplerian motion, we calculate 
the bisector velocit y spans (BVS) for spectral lines outside the I2 absorption region. As described in 
Brown et al.l (|2008l ). this calculation offers information regarding the shapes of the stellar absorption 
lines in our spectra. As photospheric activity that might influence RV measurements occurs, it alters 
the shapes of these lines, causing a corresponding shift in the BVS. Our BVS, then, offer a record 
of the activity of our target stars. Provided our planetary signals are real, we expect the BVS to 
be uncorrelated with our RV measurements. For each RV point, we have a corresponding BVS 
measurement computed from the average of the stellar lines outside the I2 region. 



^available at 



http://www.as.utexas.edu/~chris/moog.html 
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3.2. Orbit Fitting 



To determine the orbital parameters of each plan etary system, we first analyz e each RV set 
using the fully generalized Lomb-Scargle periodogram (jZechmeister Kiirsteril2009l ). We estimate 
the significance of each peak in the periodogram by assigning it a preli r ninary false-alarm prob- 
ability (FAP) according to the method described in ISturrock &: Scargld ([20101) • To conclude our 
periodogram analysis, we assign a final FAP for each planet using a boo tstrap resampl i ng te chnique. 
Our bootstrap method, which is analogous to the technique outlined in lKiirster et alj (|l997l ) retains 
the original time stamps from the RV data set, and selects a velocity from the existing set (with 
replacement) for each observation. We then run the generalized Lomb-Scargle periodogram for the 
resampled data. After 10,000 such trials for each data set, the FAP is taken to be the percentage 
of resampled periodograms t hat produced higher power than the original RV set. Power in the 
Zechmeister &: Kiirsteii (j2009l ) periodogram is given by Ax^ = Xo ~ Xpi or the improvement of fit 
of an orbit at period P over a linear fit, so a "false positive" result in the bootstrap trial occurs 
when a random sampling of our measured RVs is better suited to a Keplerian orbit than the actual 
time series. In Table El we include both FAP estimates for comparison, but we adopt this bootstrap 
calculation as our formal confidence estimate. 



The periods identified in the periodograms are then used as initial esti mates for Keplerian 
orbital fits, which we perform with the GaussFit program ([Jeffervs et al.lll988l ). We list the orbital 
parameters of each planet in Table [Gj and describe the individual sy stems in detail below . As a 
consistency check, we also fit each orbit using the SYSTEMIC console (jMeschiari et al.ll2009l ). find- 
ing good agreement between our results in every case. An additional advantage of the SYSTEMIC 
console is that it computes a stellar "jitter" term-a measure of random fluctuations in the stellar 
photosphere-(e.g. iQueloz et al.ll200ll ) for each star, which we include in Table El 



4. New Planetary Systems 
4.1. HD 79498 

4-. 1.1. RV Data and Orbit Modeling 

Our RV data for HD 79498 consist of 65 observations taken over 7 years from the 2.7 m 
telescope. The data have an RMS scatter of 18.3 m/s and a mean error of 4.15 m/s, indicating 
signiflcant Doppler motion. Our RVs are listed in Table (H and plotted in Figure [H 

The periodogram for the HD 79498 RVs reveals a large peak around 1815 days, with a pre- 
liminary FAP of 4.09 X 10^^. Performing a Keplerian fit with P = 1815 days as an initial guess, 
we find a single-planet solution with the parameters P = 1966 days, e = 0.59, and K = 26.0 m/s, 
indicating a planetary companion with Msini = 1.34Mj at a = 3.13 AU. Considering the high 
eccentricity and the width of the periodogram peak, these values are in good agreement with our 
initial guess. This fit produces a reduced = 1-77 with an RMS scatter of 5.13 m/s around the 
fit. We note that we have not included the stellar "jitter" term in our error analyses, but we have 
verified that most of the excess above unity can be attributed to our fitted value of 2.76 m/s. 
In a series of 10^ trials of our bootstrap FAP analysis, we did not find an improvement in Ax^ for 
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any resampled data set, resulting in an upper limit of 10~^ for the FAP of planet b. Fi gure [1] shows 
our fit, plotted over the time-series RV data. We include the full parameter set for HD 79498b in 
Table H 

Our analysis of HD 79498b admits a second, slightly different orbital solution of nearly equal 
significance. This solution converges at P = 2114 days, e = 0.61, and K = 26.2 m/s. The corre- 
sponding planetary parameters then become Msini = 1.35Mj and a = 3.27 AU. For comparison, 
we include the plot of this fit in Figure [TJ Although the shorter-period solution may appear to be 
driven mainly by the first data point, we note that our fitting routines converge to the 1966-day 
period regardless of whether that point is included, and that the fitting statistics are still better 
for the shorter period with the point excluded. Clearly, the qualitative properties of the planet 
remain unchanged regardless of the choice of parameters. However, because of the slightly better 
fitting statistics (x^ = 1-77 versus = 1-84, RMS = 5.13 m/s versus RMS = 5.50 m/s) for the 
first model, we adopt it as our formal solution. 

We computed the periodogram of the residuals around the fit of HD 79498b to search for 
additional signals. We see no evidence of additional planets in the system. We note that no 
additional signals appear in the residuals of the alternate fit for planet b either. 



4-1-2. Stellar Activity and Line Bisector Analysis 

With a log R'jjj^ = —4.66, HD 79498 appears to be a low-activity star, and our indicators 
corroborate that notion. Our line bisector velocity spans are well behaved, with an RMS scatter 
of 19 m/s, below the K amplitude of planet b. The BVS are uncorrelated with the measured RVs, 
and a periodogram reveals no periodicity in the bisector velocities. Likewise, the Shk time series 
shows no significant signals, and no correlation with the RV series. All measures suggest that HD 
79498 is a very quiet star, and should neither mask nor artificially produce large-amplitude RV 
signals such as the one discussed above. 



4.1.3. Stellar Companions 



HD 79498's location on the sky places it near two faint {V = 10, 11) stars, each separated by 
approximately 60 arcseconds (jDommanget &: Nysll2002l ). When examining the proper motions of all 
three stars, though, we see that only the souther n companion, Bp-|-23 2063B, is actually associated 
with HD 79498, making it a double star system (jBonnarel et al.ll2000l ). At a distance of 49 parsecs, 
the companion star is located at a minimum distance d ~ 2900 AU from HD 79498. Given an overly 
generous mass estimate of M = IMq, the secondary star would impart a gravitational acceleration 
of just GMoi/d'^ = 0.022 m s~^ yr~^, well below the sensitivity of our instrument. Furthermore, at 
2900 AU, the orbital period of the companion would be sufficiently long to appear as a linear slope 
in our data, not the orbit disc ussed above . It is possib le, however, that the presence of this distant 



star may excite Kozai cvc les (IKozai 



high eccentricity (see, e.g. iKatz et al 



1962 



2011 



Lidovlll962l) in HD 79498 b, which serve to maintain its 



Wittenmver et al.ll2007l ) 
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4.2. HD 197037 

4.2.1. RV Data and Orbit Modeling 

We have obtained 113 RV points for HD 197037 over 10 years from the 2.7 m telescope. These 
data have an RMS scatter of 13.1 m/s and a mean error of 7.65 m/s. We report our measured RVs 
in Table [3l and plot them as a time series in Figure [2j 



Our periodogram analysis of HD 197037, shown in Figure 3(a) indicates a strong peak around 
P ~ 1030 days. The power in this peak corresponds to a preliminary FAP of 8.53 x 10^^. Our 
one-planet Keplerian model yields parameters of P = 1036 days, e = 0.22, and K = 15.5 m/s, 
showing excellent agreement with the prediction of our periodogram. The inferred planet has a 
minimum mass of 0.79Mj and lies at an orbital separation of 2.07 AU. We overplot our model with 
the RVs in Figure [2l 

The addition of a single Keplerian orbit only reduces the residual RMS of our RVs to 9.18 m/s, 
which is nearly a factor of two worse than our other one-planet fits. Furthermore, the periodogram 



of the residual RVs (Figure 3(a) ) shows a significant increase in power at long periods. While any 
additional periods are too far outside our observational time baseline to properly evaluate, we can 
generate a preliminary two-planet fit with a second ~ 0.7Mj planet with a period around 4400 days 
and an eccentricity 0.42. However, our current RV set can also be modeled as a single planet and a 
linear trend with slope —1.87 ± 0.3 m s~^ yr~^. Both fits give a of 1.10 and an RMS scatter of 
8.00 m/s, so we adopt the more conservative planet-plus-slope model pending further observations. 
HD 197037 has no known common proper motion companions within 30 arcseconds, so it is very 
possible that we are seeing evidence for a distant giant planet or brown dwarf companion. The final 
orbital parameters reported in Table [6] are derived from the planet+slope model. We note that we 
see no additional signals in the residuals of either the two-planet or planet-slope fits. 



4.2.2. Stellar Activity and Line Bisector Analysis 

Even after accounting for the linear trend in the residuals around our fit to planet b, the scatter 
in our RVs is still higher than we typically expect from our 2.7 m data. Because HD 197037 is an 
earlier type (F7) than the other stars discussed in this paper, we expect a lower precision as a result 
of fewer spectral lines to determine velocities. Additionally, with logi?^j^ = —4.53, it is the most 
active of the stars presented here. With an RMS scatter in the bisector velocity spans of 13 m/s, we 
must take particular care to ensure we have not mistaken an activity cycle for a planetary signal. 
In Figure m we show our measured BVS and Shk indices, which we use to evaluate the influence 
of stellar photospheric activity on the signal of HD 197037b. Using both a Pearson correlation test 
and least-squares fitting, we find no correlation to the RVs for the BVS. While the Shk indices 
appear somewhat correlated with the RVs to the unaided eye, the Pearson correlation coefficient 
of -0.26 indicates that the relation is not statistically significant. Periodogram analysis (Figure 



3(b) ) shows no power at the 1030 day peak in the BVS or Shk time series. We do see a modestly 
(FAP = 0.01) significant peak at 19.1 days in the periodogram of our Shk measurements, which we 
speculate may be the stellar rotation period. Although we are admittedly unable to completely rule 
out the possibility of stellar activity as the source of our observed RV signal, the lack of correlation 
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between our activity indicators and the velocities, and the absence of periodicity in Shk and BVS 
around the fitted period lead us to the conclusion that a planetary orbit is the most likely cause. 

While the stellar activity measurements reinforce the planetary nature of the primary RV 
signal of planet b, the analysis of the long-term trend is not so clear. Our periodograms of the 
BVS and Shk series both show an increase in power at very long periods, matching the behavior 
of the trend in the residual RVs, albeit at much lower power. Furthermore, we see a correlation 
between Shk and the residual RVS around the one-planet fit with a Pearson correlation coefficient 
of -0.33, which is significant for the size of our data set. The residual RVs are uncorrelated with 
the BVS, though, so the trends may be coincidental. HD 197037 will need continued monitoring 
to determine the true nature of this long-period signal. 

4.3. HD 220773 

4.3.1. RV Data and Orbit Modeling 

Our RV data set (Table HD for HD 220773 consists of 43 HET/HRS spectra taken over 9 years 
between July 2002 and August 2011. The data have an RMS of 11.7 m/s and a mean error of 4.11 
m/s. 

Based on our data, we find evidence for a highly eccentric giant planet on a long-period orbit. 
Because of its high eccentricity, the period of planet b does not appear at significant power in our 
periodogram analysis. However, the high RMS of our RV set, combined with the characteristic 
eccentric turnaround of the velocity time series (Figure [5]), strongly indicate the presence of a 
planetary companion. In order to offset the lack of information from the periodogram, we have 
run GAUSSFIT with a broad range of periods (2500d-4500d) and eccentricities (0.4-0.7). All fits 
converge unambiguously to a period of 3725 days with an eccentricity of 0.51, indicating a 1.45-/Vfj 
planet at 4.94 AU. We include the plot of this model in Figure [5l The final fit gives a reduced 
of 3.14 and an RMS scatter of 6.57 m/s. Because the reduced is so high, we have computed the 
model again after adding the 5.10 m/s "jitter" term in quadrature to the errors listed in Tabled 
While our fitted parameters and uncertainties do not change significantly from the values given in 
Table m the reduced drops to 1.18, lending additional confidence to our solution. We do not 
currently see any evidence for additional signals in the residual RVs. 

To confirm the orbital parameters listed in Table [6l we use a genetic algorithm to explore 
the parameter space and evaluate the likelihood of the null hypothesis. The algorithm fits a grid 
of parameters (number of planets, masses, periods, eccentricities) to the RV set, allowing for a 
thorough exploration of how behaves in response. Areas of parameter space which do not match 
the data are iteratively rejected, allowing the routine to converge on an optimal solution. We 
have performed 10,000 iterations of the algorithm with our RV data, considering periods between 
3000 and 10,000 days. For this experiment, we have again added the stellar "jitter" term to our 
measurement errors. As shown in Figure El the genetic algorithm reaches x^ < 1-20 for a one-planet 
model with periods close to our fitted value in Table [H This is a dramatic improvement over a 
zero-planet model, which yields = 3.65. 
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4-3.2. Stellar Activity and Line Bisector Analysis 

Our bisector velocity spans for HD 220773, which have an RMS scatter of 24 m/s, show no 
correlation to our measured RVs, and a periodogram analysis of the BVS shows no power around 
the period of planet b. While the periodogram signal for the RVs did not meet the criterion for a 
positive detection, there was a broad peak centered around the ~ 4000 day period of the planet. 
The BVS, on the other hand, show no evidence whatsoever for a long-period trend or signal. We 
conclude, then, that the BVS do not indicate stellar activity which could mimic the behavior of this 
long-period planet. Likewise, while the HRS does not provide time-series Shk information, our 
2.7 m stellar spectrum does offer a "snapshot" value of logi?^^. From this spectrum, we derive 
logi?'^^ = —4.98, which is consistent with the low activity level suggested by the BVS analysis. 



5. Updated Planetary Parameters for HD 155358 



HD 155358 has pre viously been identifie d as a two-planet system by the McDonald Observatory 
Planet Search Program (jCochran et al.ll2007l ) . It is a notable system because it hosts two giant plan- 
ets despite having a measured [Fe/H] among the lowest of any stars with subste llar companions of 
an unambiguously planetary na ture (excepting HIP 130 4 4 (Setiawan et aklboill ). which apparently 
originates outside the Galaxy). iFuhrmann &: Bernkopfl (120081 ) analyze its chemical abundances in 
more detail, leading them to claim that the star is actually a mem ber of the thick disk popu lation of 
the Galaxy, a claim strengthened by the > 10 Gyr age estimate of ICasagrande et al.l (|201ll ). While 
our updated analysis suggests a metallicity of [Fe/H] = -0.51, rather than our original estimate of 
[Fe/H] = -0.68, the system is still extremely metal-poor relative to the other known exoplanetary 
systems. Continued monitoring of HD 155358 has caused us to reevaluate our prev iously reported 



orbita l parameters, resolving the ambiguity in the period of planet c discussed in ICochran et al 
(I2OO7I ). 



5.1. RV Data and Orbit Modeling 



Since the publication of ICochran et al.l (j2007l ). we have obtained 51 additional RV points for 
HD 155358 from HET. We have re-reduced all of our spectra with the latest version of AUSTRAL 
for consistency, and include all 122 points in Table [2l The entire set of velocities, shown as a time 



series in Figure 7(a) , covers approximately 10 years from June 2001 until August 2011. These RVs 



have an RMS scatter of 30.2 m/s, much higher than we expect from a mean error of just 5.38 m/s. 



We include the periodogram and window function for HD 155358 in Figure 
we see a significant peak around 196 days, with a preliminary FAP of 1.15 x 10 
one-planet orbit from this peak, we find an RV amplitude K = 32 m/s with a period P 
and an eccentricity of e = 0.23. Based on this fit, HD 155358b has a minimum mass of 0.85 Jupiter 
masses at a = 0.64 AU. 



As expected. 
Fitting a 
195 days 



Cochran et al.l (|2007l ) note that the periodogram for the residuals around the one-planet fit to 



planet b initially showed power around 530 days and 330 days. Comparing fits at both periods 
produced better results at the longer period for our data at that time. However, the periodogram of 
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the residuals around planet b for our current data (Figure 8(a) ) clearly indicates the shorter period 



as the true signal. The peak at 391 days has a preliminary FAP too low for the precision of our 
code (approximately 10~^^), and we no longer see additional peaks at longer periods. While the 
window function for our sampling does show some power at the one-year alias, the period of this 
periodogram peak is sufficiently separated from the yearly alias (Figure [8(b)[ ) that we are confident 
the signal is not due to our sampling. 

Using 391 days as a preliminary guess for the period of planet c, we performed a two-planet fit to 
our RVs. This updated fit changes the parameters of planet b to P = 194.3 days, Msini = 0.85Mj, 
and e = 0.17, with the orbital separation remaining at 0.64 AU. Planet c then converges to a period 
of 391.9 days, with e = 0.16 and K = 25 m/s. The derived properties for planet c then become 
M = 0.82Mj and a = 1.02 AU. We plot this orbit over our RV data in Figure |7(a)] The addition 
of planet c reduces our RMS to 6.14 m/s with a reduced of 1.41. As a consistency chec k, we 



test a fit with parameters for planet c more closely matching those from lCochran et al.l ( 20071 ) . but 
find no satisfactory solution at this longer period. Holding the parameters for planet c fixed at the 
values derived in that study result in a reduced of 11.2 and a residual RMS of 15.4 m/s. 



5.2. Stellar Activity and Line Bisector Analysis 

The BVS for HD 155358 display an RMS scatter of 25 m/s, and are uncorrelated with both 
the RVs and the residuals to the 2-planet fit. Additionally, we see no significant peaks in the 
periodogram of the BVS time series. Also, while we present HET velocities here, we do have some 
2.7 m spectra from a preliminary investigation of HD 155358, allowing us to examine the Ca H and 
K indices. Shk shows very little activity, and we derive a logi?^^ of -4.54. The hypothesis that 
the two large-amplitude signals observed here are due to stellar activity can therefore conclusively 
be ruled out. 



5.3. Dynamical Analysis of 3-body System 



Even at the separation claimed in ICochran et al.l (|2007l ). planets b and c are close enough to 
interact gravitationally. Our updated orbital model now indicates they are actually much closer 
together. While their periods suggest the planets are in a 2:1 mean motion resonance, our prelimi- 
nary orbital simulations showed the system to be unstable for a range of input values. We therefore 
decided to perform a highly detailed dynamical study of the system to investigate whether the 
orbits that best fit the data are indeed dynamically feasible. To do this, we performed over 100,000 
unique simulations of the H D 155358 systern . using the Hybrid integrator within the n-body dy- 
namics package MERCURY (Ichamberslll999l ). 



To systematically address th e stability of t he HD 155 358 system as a funct ion of the orbits 
of planets b and c, we followed iHorner et al.l (|201ll ) and [Marshall et al.l (|2010l ). and examined 
test systems in which the initial orbit of the planet with the most tightly constrained orbital 
parameters (in this case planet b) was held fixed at the nominal best fit values. The initial orbit of 
the outermost planet was then systematically changed from one simulation to the next, such that 
scenarios were tested for orbits spanning the full ±3(7 error ranges in semi-major axis, eccentricity. 
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longitude of periastron and mean anomaly. Such tests have already proven critical in confirming 



or rejecting planets suggested to move on unusual orbits (e.g. IWittenmver et al.ll2011bl ). and allow 



the construction of detailed dynamical maps for the system in orbital element phase space. 

Keeping the initial orbit of the innermost planet fixed, we examined 31 unique values of semi- 
major axis for planet c, ranging from 0.96 AU to 1.08 AU, inclusive, in even steps. For each 
of these 31 initial semi-major axes, we studied 31 values of orbital eccentricity, ranging from the 
smallest value possible (0.0) to a maximum of 0.46 (corresponding to the best fit value, 0.16, plus 
three sigma). For each of the 961 a-e pairs, we considered 11 values of initial mean anomaly and 
initial longitude of periastron (w), resulting in a total suite of 116281 (31 x 31 x 11 x 11) plausible 
architectures for the HD 155358 system. 

In each simulation, the planet masses were set to their minimum (M sin i) values; the mass of 
the innermost planet was set to 0.85Mj, and that of the outermost was set to 0.82Afj. As such, 
the dynamical stability maps obtained show the maximum stability possible (since increasing the 
masses of the planets would clearly increase the speed at which the system would destabilize for any 
nominally unstable architecture). The dynamical evolution of the two planets was then followed 
for a period of 100 million years, or until one of the planets either collided with the central star, 
was transferred to an orbit that took it to a distance of at least 10 AU from the central star, or 
collided with the other planet. 

The results of our simulations are shown in Figures [911101 We present the multi-dimensional 
grid of orbital parameters over which we ran our simulations as two-dimensional cross sections, 
indicating the mean and median stability lifetimes over each of the runs at each grid point. So, 
for example, each grid point in our plots in a-e space reveals either the mean or the median of 121 
unique simulations of an orbit with those particular a-e (i.e. 11 in times 11 in mean-anomaly). 
Similarly, each point in the a-co plots corresponds to the mean (or median) of 341 separate trials 
(31 in e times 11 in uj). 



The first thing that is apparent from Figure 9(a) , which shows the mean lifetime of the system 



as a function of semi- major axis and eccentricity, is that the stability of the system is a strong 
function of eccentricity. Solutions with low orbital eccentricity are typically far more dynamically 
stable than those with high eccentricity. In addition, the stabilizing effect of the 2:1 mean-motion 
resonance between planets b and c can be clearly seen as offering a region of some stability to even 
high eccentricities between a ~ 1.05 and 1.25 AU. However, it is apparent when one examines the 
median lifetime plot (Figure |9(b)[ ) that a significant fraction of eccentric orbits in that region are 
dynamically unstable (the reason for the apparently low median lifetimes in that region) . The reason 
for this is that the dynamical stability of orbits in this region, particularly for high eccentricities, is 
also a strong function of the initial longitude of periastron for planet c's orbit. Indeed, comparison 
of this figure to those showing the influence of the longitude of periastron for planet c reveals 
that the most stable regions therein lie beyond the la error bars for the nominal orbit. For low 
eccentricities (e < 0.1), in the vicinity of the 2:1 resonance, then the orbit is stable regardless of the 
initial longitude of periastron, but for high eccentricities (e > 0.15), the stable regions lie towards 
the edge of the allowed parameter space, making such a solution seem relatively improbable. 



It is apparent from close examination of Figure |9(b) (which shows the median stability as a 



function of semi- major axis and eccentricity) that the key determinant of the stability of the system 
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(particularly for non-resonant orbits) is actually the periastron distance of planet c. The sculpted 
shape of the stability plot, outwards of a ~ 0.98 AU, i s very similar to that observed for the propose d 



planets around the cataclysmic variable HU Aquarii (jHorner et al.ll201ll : IWittenmver et al.ll2011bl ). 
As was found in that work, the dividing line between unstable and stable orbits again seems to fall 
approximately five Hill radii beyond the orbit of the innermost planet. Any orbit for planet c that 
approaches the orbit of b more closely than this distance will be unstable on astronomically short 
timescales (aside from those protected from close encounters by the effects of the 2:1 mean motion 
resonance between the planets). On the other hand, orbits which keep the two planets sufficiently 
far apart tend to be dynamically stable. We note, too, that the sharp inner cut-off to this broad 
region of stability, at around 0.99 AU, corresponds to the apastron distance of the innermost planet 
plus five Hill radii. Once again, orbital solutions that allow the two planets to approach closer than 
five times the Hill radius of the inner planet destabilize on relatively short timescales. 

The results of our dynamical simulations show that a large range of dynamically stable solutions 
exist within the la errors on the orbit of HD 155358c. Given the breadth of the resonant feature 



apparent in Figure 9(a) (the mean lifetime as a function of a and e), it seems clear that all orbits 
within la of the nominal best fit will be strongly influenced by that broad resonance. As such, 
it seems fair to conclude that these two planets are most likely trapped within their mutual 2:1 
mean-motion resonance, although it is also apparent that non-resonant solutions also exist that 
satisfy both the dynamical and observational constraints. The stability restrictions are consistent 
with our orbital fit, but the parameters listed in Table [6] have not been modified to include the 
information obtained from the orbital simulations. 



5.4. Habitability of Exomoons 



Using the [Fe/H l and Tpjf derived herei n, an [a/Fe] of 0.32 (jFuhrmann &: Bernkopal2008ll . and 



an ae 



ge of 10.7 Gyr (jCasagrande et al 



2011 



, we have fit Yonsei- Yale isochrones (IDemaraue et al 



20041) to HD 155358. Wi th mass estimates ranging from O.87M0 (jLambert et al. 



1991 



) to O.92M0 
1.14L0 and 



(jCasagrande et al.l 1201 ll ). HD 155358 has a luminosity somewhere between L^, 

= I.GTLq. If we assume the location of the habitable zone scales as \/L^, then at an or- 
bit al separation of 1.06 AU, planet c lies within the habitable zone of its parent star (as defined 
by iKasting et al.l Il993l). Typically, Jo vian planets are not considered to be potential habitats 
for Earth- like life (|Lammer et al.ll2009l ). However, it has been suggested that such planets could 



host potentially habitable satellites (e.g. iPorter Grundvll201ll ). or even that they could dynami- 
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cally capture plane t-sized objects during their migration as satellites or Trojan companionj^ (e.g. 



Tinnev et al. 



2011 



). In fact, the presence of two giant planets on relatively close-in orbits may 
indirectly in crease the water cont ent on terrestrial satellites through radial mixing of planetesimals 
rich in ices (IMandell et al.l 120071 ) . a key requirement for such objects to be considered habitable 
(e.g. iHorner &: Jonesll2010l ). Furthermore, any satellites sufficiently large to be considered habitable 
would also be subject to significant tidal heating from their host planet, which would likely act to 
increase their habitability when they lie towards the outer edge of Kasting's habitable zone (as they 
would have earlier in the life of the syst em). The induced tect onics would also potentially improve 



the habitability of any such moons (e.g. iHorner &: Joneal2010l ). 



Our dynamical analysis suggests that the two planets in the HD 155358 system are most likely 
trapped in mutual resonance. It might initially seem that such orbits would be variable on the long 
term, such that they would experience sufficiently large excursions as to render any satellites or 
Trojan companions uninhabitable. However, we know from our own Solar system that long term 
resonant captures can be maintained on timesca les comparable to the age of the Solar system (e.g. 
the Neptune Trojans. iLykawka &: Horned l2010l ). As such, it is reasonable to assume that, if the 
planets are truly trapped in mutual resonance, they could have been on their current orbits for at 
least enough time for any moons or Trojans of planet c to be considered habitable. However, the 
current best-fit eccentricity for planet c is suffi ciently high that it l ies at the outer limit of the range 
in e that would allow for a habitable exomoon ([Tinnev et al.ll201ll ). which might limit the potential 
habitability of any exomoons in the system. That said, the one-sigma range of allowed orbits 
extends to relatively low eccentricities (e = 0.06), which is certainly compatible with potential 
habitability. As such, it is possible that such moons, if present, could be potential habitats for life. 



6. Discussion 

Of the planets discovered via the radial velocity method, only around 7 percent have a > 3 
AU. The addition of the planets presented here therefore represent a significant contribution to 
that sample. As more such objects start to fall within the detection limits of RV surveys, our 
findings provide interesting comparison cases to begin to look for trends in the long-period gas 
giant population. 

Continued monitoring of the HD 79498, HD 197037, and HD 220773 systems is crucial for 



'^Planetary Trojans are a particularly fascinating population of objects, trapped in 1:1 mean-motion resonance 
with their host planet. Typically, stable Trojans follow horseshoe-shaped paths that librate around either the 
leading- or trailing-Lagrange points in that planets orbit, located sixty degrees ahead a nd behind the planet (for a 
good illustration of such orbits, we direct the interested reader to Horner fc Evans 20061 ). 



Within our Solar system, both J upiter and Neptune host significant popu lations of Trojans that were cap- 

For those planets, the captured 



tured during their migration (e.g. iMorbidelli et al 



2005 



Lvkawka et al 



200' 



i). 



Trojans are typically small, but there is nothing to prevent a giant planet from capturing an Earth-mass object as a 
Trojan during its migration. Once c aptured, and once the migration stops, such ob jects can be dynamically stable 
on timescales of billions of years (e.g. Horner fc Lvkawka 20ld : Lvkawka et al. 2011), even when moving on orbits of 
significantly different eccentri city and inclination to their host planet. Whilst the detection of such planets would 
no doubt be challenging (e.g. iFord fc HolmanI 120071 ). they remain an intriguing option in the search for habitable 
worlds. 
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refinement of planet formation and migration theories. The presence or absence of ad ditional smaller 



planets will shape om' understanding of the migrational history of these systems. iMandell et al 



(j2007l ). for example, predict markedly different outcomes for the formation and water content of 
sub-giant planets after migration of a Jovian planet depending on whether or not gas drag plays a 
significant role in planetary migration. 

Of the four stars presented herein, at least one-HD 155358-hosts multiple gas giant planets. 
Additionally, HD 197037 shows tentative evidence of a yet-undetected substellar companion at 
a ~ 5.5 AU. If the occurrence of multiple- Ju piter systern s is as common as it appears in this very 
small sample, it would support the claim bv ISumi et al.l (j201ll ) that the large population of free- 
floating planetary-mass objects within the Galaxy form i n protoplanetary disks as planets. On the 
other hand, current RV results (IWittenmyer et al.ll2009l ) indicate that multi-giant planet systems 
should actually be quite rare, at least inside of 2-3 AU. Understanding whether that result remains 
valid for planets further out, though, will rely heavily on detections of residual long-term trends 
for planets hke HD 197037b. 

The orbital evolution of the HD 155358 system is of particular interest, as it provides a compar- 
ison case for theoretical w ork exploring the formation and migration of the Jupiter-Saturn system 
(e.g. iTsiganis et al.l l2005l ) . While Jupiter and Saturn appear to have crossed the 2:1 MMR and 
later separated, the HD 155358 giants remain locked in resonance, presumably for the entire ~ 10 
Gyr lifetime of the system. Accounting for the difference in these final configurations most likely 
requires different initial architectures, disk masses, and encounter histories for the two systems. 

While the stability simulations presented here offer some constraints on the orbital configu- 
rations of the HD 155358 system, the true geometries of those planets' orbits will only be fully 
understood if their mutual inclinations are measured. Unfortunately, the HD 155358 planets will 
only be accessible to the next generation of astrometric instrumentation, as their predicted as- 
trometric displace ments are below the 0.2 mas precision limit of the HST Fine Guidance Sensor 
(jNelan et al.ll20ld ). On the other hand, HD 220773b, which has a predicted displacement of 0.242 
mas, might be an interesting astrometry target for FGS, both for the purpose of determining a true 
mass, and to search for outer companions to this distant planet. 

The refined orbital parameters for HD 155358c places it, to zeroth order, within the habitable 
zone. While the planet itself is likely inhospitable to any Earthlike lifeforms, the possibility of 
habitable m oons or Trojans make it an i nteresting datum for examining more exotic environments 



for biology (jSchulze-Makuch et al.l 120111 ). Given the advanced age of the system, any intelligent 
life residing on planet c is potentially far more advanced than our own civilization. The techno- 
logical advances accompanying such extended development may make interstellar broadcasts or 
beac ons energetically and f inancially feasible, mak ing HD 155358 an interesting target for SETI 
(e.g. I J. Benford et ahlboiol : Ig. Benford et ahlboiol ) 



Our results are somewhat atypic al when viewed iii the c ontext of the metallicity-frequency 



correlation for giant planets (see, e.g. iFischer &: Valentil l2005l ). Only HD 79498 has significantly 
super-solar metallicity, while HD 197037 and HD 155358 qualify as metal-poor. In the case of 
HD 197037, its higher mass and earlier spectral type increa se its likelihood of forming Jovian 
planets, at least somewhat offsetting any metallicity effects (jjohnson et al.ll201ll ). HD 155358, 



on the other hand, remains a true anomaly despite the slightly higher [Fe/H] reported herein. 
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While the sample presented here is obviously too small to make even tentative statements about 
the validity of observed correlations, it will be interesting to see whether gas giants continue to be 
found preferentially around metal-rich stars as long-term RV surveys begin to reveal a large number 
of Jupiter-mass planets. If the removal of the period bias on the giant planet census reveals a large 
population of Jovian planets around metal-poor stars, it will serve as strong evidence that at least 
some gas giants form through gravitational instability of the protoplanetary disk rather than core 
accretion (lBosal2002l ). Furthermore, the discovery of thick disk planets suggests a planet formation 
history essentially spanning the age of the Galaxy. 



The authors wish to thank Anita Cochran, Candace Gray, and Diane Paulson for their contribu- 
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2453500 2454000 2454500 2455000 2455500 2456000 
Barycentric Julian Date (days) 

Fig. 1. — Top: Radial velocity data for HD 79498. The best-fit orbit model is shown as a solid 
red line. Our second acceptable fit is also included, shown as a dashed blue line. Middle: Residual 
RVs after subtraction of a one-planet fit. Bottom: Shk as measured at each RV point. 
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2452000 2453000 2454000 2455000 2456000 

Barycentric Julian Date (days) 

Fig. 2. — Top: Radial velocity data for HD 197037. The best-fit orbit model is shown as a solid 
red line. Middle: Residual RVs after subtraction of a one-planet fit and a linear trend. Bottom: 
Shk as measured at each RV point. 
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Fig. 3. — a. Prom top: [1] Generalized Lomb-Scargle periodogram for HD 197037 RVs. [2] The 
same periodogram for the residual RVs after subtracting a one-planet fit. [3] Periodogram of the 
residual RVs after subtracting a one-planet fit and a linear trend. [4] Periodogram of our time 
sampling (the window function). 



b. Generalized Lomb-Scargle periodograms for the BVS {top) and Shk indices (bottom) of 
our spectra for HD 197037. The dashed lines indicate the approximate power level for a FAP of 
0.01. 
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Fig. 4. — Left: Bisector velocity spans plotted against our measured RVs (top) and residual RVs 

to a one-planet fit (bottom) for HD 197037. Right: Shk indices plotted against our measured 
RVs {top) and residual RVs to a one-planet fit (bottom) for HD 197037. The error bars shown are 
representative of the data set. 
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Fig. 5. — Top: Radial velocity data for HD 220773. The best-fit orbit model is shown as a red line. 
Bottom: Residuals to a one-planet fit. 
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Fig. 6. — Genetic algorithm results for HD 220773. Each circle indicates the minimum reached 
by fitting a one-planet orbit at the given period. 
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b. Phase plots for planets b (top) and c {bottom). In each plot, the signal of the planet 
not shown has been subtracted from the RVs. 
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Fig. 8. — a. From top: [1] Generalized Lomb-Scargle periodogram for HD 155358 RVs. [2] The 
same periodogram for the residual RVs after subtracting a one-planet fit. [5] Periodogram of the 
residual RVs after subtracting a two-planet fit. [4] Periodogram of our time sampling (the window 
function). The dashed lines indicate the approximate power level for a FAP of 0.01. 



b. Periodogram for planet c (black) with the window function (red). Note that the two 
peaks do not overlap. 
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Fig. 9. — Results of our stability analysis for the HD 155358 planetary system in a — e space. At 
each grid point, we have run 25 dynamic simulations with planet c at that point, with planet b 
remaining fixed at the configuration given in Table [6l The resulting intensity indicates the mean 
(a) or median (b) lifetime of the planet over those 25 runs in log years. Longer lifetimes indicate a 
more stable configuration. The crosshairs indicate our best-fit parameters as derived from the RV 
data (Table ED. 
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Fig. 10. — Results of our stability analysis for the HD 155358 planetary system in a — a; space. At 
each grid point, we have run 25 dynamic simulations with planet c at that point, with planet b 
remaining fixed at the configuration given in Table [6l The resulting intensity indicates the mean 
(a) or median (b) lifetime of the planet over those 25 runs in log years. Longer lifetimes indicate a 
more stable configuration. The crosshairs indicate our best-fit parameters as derived from the RV 
data (Table ED. 
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Table 1: Radial Velocities for HD 79498 
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TableUl cont'd. 


BJD - 2450000 


Radial Velocity 


Uncertainty 


Shk 




(m/s) 


(m/s) 




5278.74248800 
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Table 2: Radial Velocities for HD 155358 



BJD - 2450000 


Radial Velocity 


Uncertainty 




(m/s) 


(m/s) 


2071.90483700 


15.90 


5.14 


2075.88682700 
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Table\M cont'd. 


BJD - 2450000 


Radial Velocity 


Uncertainty 




(m/s) 


(m/s) 
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3.17 
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5.67 
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Table\M cont'd. 


BJD - 2450000 


Radial Velocity 


Uncertainty 




(m/s) 


(m/s) 


4369.61238300 


8.15 
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3.28 
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Table 3: Radial Velocities for HD 197037 
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Table\^ cont'd. 


BJD - 2450000 


Radial Velocity 


Uncertainty 


Shk 




(m/s) 


(m/s) 




4346.72043545 


0.46 


7.49 


0.1655 ±0.0229 


4375.64708486 


-11.15 


7.87 


0.1714 ±0.0230 


4375.65745814 


-10.85 
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0.2151 ±0.0265 
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0.1685 ±0.0235 
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0.1528 ±0.0223 
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-3.48 


7.97 


0.1618 ±0.0214 


4405.56534582 


-10.03 


7.13 


0.1655 ±0.0231 


4570.97352687 


19.43 


8.12 


0.1701 ±0.0243 


4604.92111685 


8.57 


8.82 


0.1820 ±0.0264 


4605.91057844 


12.66 


7.48 


0.1732 ±0.0255 


4662.95425192 


11.46 


8.31 


0.1771 ±0.0251 


4662.96610624 


-1.38 


8.21 


0.1773 ±0.0251 


4663.81542937 


19.18 


9.20 


0.1851 ±0.0256 


4663.85662433 


2.50 


7.23 


0.1725 ±0.0244 


4663.88867894 


4.10 


6.74 


0.1691 ±0.0235 


4663.90391404 


10.16 


6.00 


0.1717 ±0.0252 


4665.78413890 


11.48 


7.78 


0.1816 ±0.0257 


4665.79430438 


-3.10 


9.86 


0.1812 ±0.0252 


4665.80583621 


0.63 


9.02 


0.1786 ±0.0242 


4731.66600573 


1.60 


11.50 


0.1734 ±0.0230 


4731.67754085 


0.72 


7.35 


0.1657 ±0.0226 


4731.68907100 


-5.82 


8.82 


0.1795 ±0.0236 


4732.75396109 


2.27 


7.88 


0.1682 ±0.0242 


4732.76548879 


0.12 


7.60 


0.1667 ±0.0223 


4732.77701662 


-2.03 


6.90 


0.1674 ±0.0238 


4733.69703238 


9.44 


8.14 


0.1701 ±0.0242 


4733.70856019 


-4.67 


7.82 


0.1677 ±0.0235 


4733.71697119 


-7.09 


9.14 


0.1727 ±0.0236 


4747.77353788 


4.25 


5.98 


0.1746 ±0.0250 


4781.65688110 


7.76 


6.41 


0.1652 ±0.0206 


4816.61232656 


10.29 


8.72 


0.1755 ±0.0209 


4986.92370887 


-4.73 


14.24 


0.1887 ±0.0262 


5025.89140933 


-37.17 


8.23 


0.1816 ±0.0263 


5048.87005530 


-37.84 


9.19 


0.1823 ±0.0239 


5048.87535157 


-12.18 


9.07 


0.1807 ±0.0241 


5075.87456767 


-15.99 


7.03 


0.1695 ±0.0230 


5075.88610277 


-13.42 


8.76 


0.1692 ±0.0228 


5100.78063288 


-21.91 


7.32 


0.1673 ±0.0212 


5135.72812260 


-27.54 


8.80 


0.1659 ±0.0226 
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Table\^ cont'd. 


BJD - 2450000 


Radial Velocity 


Uncertainty 


Shk 




(m/s) 


(m/s) 






5172.58815669 


-4.45 


8.31 





1728 ±0.0219 


5400.78427274 


-10.86 


6.31 





1706 ± 0.0243 


5435.78861255 


-15.49 


8.10 





1714 ± 0.0246 


5469.78414305 


-13.65 


7.30 





1785 ± 0.0224 


5492.76854834 


-11.30 


7.62 





1627 ± 0.0222 


5523.68425857 


2.79 


7.30 





1601 ± 0.0208 


5527.64465815 


6.66 


8.70 





1742 ± 0.0234 


5527.65392663 


3.10 


5.77 





1669 ± 0.0223 


5547.59458469 


9.04 


8.11 





1649 ± 0.0249 


5644.98896924 


4.16 


6.79 





1519 ±0.0239 


5722.92872656 


7.82 


5.48 





1523 ± 0.0216 


5758.92483333 


6.08 


6.53 





1721 ± 0.0260 


5789.68097414 


-0.49 


6.97 





1705 ± 0.0285 


5791.86463193 


13.74 


8.20 





1768 ± 0.0263 


5792.79925378 


-0.60 


6.61 





1742 ± 0.0253 


5841.78442437 


-8.29 


7.41 





1665 ±0.0251 


5842.79215631 


-10.22 


6.66 





1649 ±0.0281 



Table 4: Radial Velocities for HD 220773 



BJD - 2450000 


Radial Velocity 


Uncertainty 




(m/s) 


(m/s) 


2479.859725 


-0.4171 


4.08 


2480.838237 


-13.5865 


4.09 


2486.820898 


-11.3423 


4.45 


2487.830499 


-16.1442 


4.38 


2948.718426 


-10.1010 


3.87 


3185.902803 


-24.8888 


5.50 


3267.686065 


-22.1946 


4.73 


3542.934309 


-35.1965 


4.76 


3604.768270 


-26.1322 


4.20 


3901.952807 


3.7467 


3.70 


3935.863371 


-5.9660 


4.19 


3979.907141 


-12.5221 


4.13 


4007.670760 


1.0665 


3.68 


4014.657584 


-0.4731 


3.14 


4015.639336 


5.8017 


4.02 


4032.599555 


14.5027 


4.60 


4044.722222 


14.3386 


4.10 


4401.744726 


12.0180 


4.51 


4428.675456 


19.4405 


3.78 


4475.548221 


14.3677 


4.52 


5114.635072 


-11.1314 


4.15 


5140.721408 


7.8337 


3.75 


5358.945550 


4.7166 


4.79 


5358.947973 


8.1944 


4.58 


5382.886175 


-5.2398 


3.81 


5382.888643 


-3.2236 


4.03 


5414.822196 


-1.2009 


3.47 
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ra6k[2 cont'd. 


BJD - 2450000 


Radial Velocity 


Uncertainty 




(m/s) 


(m/s) 


5414.824720 


-10.5412 


3.58 


5441.730451 


-0.5391 


3.39 


5441.732880 


-4.1901 


3.75 


5466.827877 


-6.5113 


4.33 


5469.822406 


-1.8755 


3.72 


5469.824889 


-1.5983 


3.69 


5532.649379 


-3.1531 


4.21 


5532.651798 


8.9105 


4.23 


5547.609328 


8.9497 


3.68 


5547.611809 


2.8300 


4.20 


5724.945915 


9.2505 


3.56 


5724.948402 


-5.5916 


3.73 


5754.866046 


-12.5525 


4.50 


5754.868528 


0.3110 


5.85 


5790.949560 


-2.0862 


3.38 


5790.952042 


-3.9868 


3.78 



Star 


Spectral 




B-V 


Mv 


Parallax'^ 


Dist. 




logs 


[Fe/H] 




Mass"* 


Age'^ 


Shk 


log R'hk 




Type 








(mas) 


(pc) 


(K) 






(km/s) 


(Mo) 


(Gyr) 






HD 79498 


G5 


8.027 


0.703 


4.58 


20.43 


49 


5740 ± 100 


4.37 ±0.12 


0.24 ±0.06 


1.27 ±0.15 


1.06 


2.70 


0.167 


-4.66 


HD 155358 


GO 


7.270 


0.550 


4.09 


23.07 


43 


5900 ± 100 


4.16 ±0.12 


-0.51 ±0.06 


0.50 ±0.15 


0.92 


10.70 


0.169 


-4.54 


HD 197037 


F7 V 


6.809 


0.519 


4.23 


30.50 


33 


6150 ± 100 


4.40 ±0.12 


-0.20 ±0.06 


1.16 ±0.15 


1.11 


1.90 


0.167 


-4.53 


HD 220773 


F9 


7.093 


0.632 


3.66 


20.60 


49 


5940 ± 100 


4.24 ±0.12 


0.09 ±0.06 


1.50 ±0.15 


1.16 


4.40 


0.159 


-4.98 



Table 5: Stellar properties. 



teharchenko Roeserl ^200% 



Casagrande et all (|201ll ) 





HD 79498b 


HD 155358 
b c 


HD 197037b'' 


HD 220773b 


Period P (days) 


1966.1 ±41 


194.3 ±0.3 


391.9 ±1 


1035.7 ± 13 


3724.7 ± 463 


Periastron Passage To 


3210.9 ± 39 


1224.8 ± 6 


5345.4 ± 28 


1353.1 ± 86 


3866.4 ± 95 


(BJD - 2 450 000) 












RV Amplitude K (m/s) 


26.0 ± 1 


32.0 ± 2 


24.9 ± 1 


15.5 ± 1 


20.0 ±3 


Mean Anomaly Mo^ 


329° ± 11° 


129° ±0.7° 


233° ± 0.9° 


186° ± 3° 


226° ± 13° 


Eccentricity e 


0.59 ± 0.02 


0.17 ±0.03 


0.16 ±0.1 


0.22 ± 0.07 


0.51 ±0.1 


Longitude of Periastron uj 


221° ± 6° 


143° ± 11° 


180° ± 26° 


298° ± 26° 


259° ± 15° 


Semimajor Axis a (AU) 


3.13 ±0.08 


0.64 ± 0.01 


1.02 ± 0.02 


2.07 ±0.05 


4.94 ± 0.2 


Minimum Mass Msini {Mj) 


1.34 ±0.07 


0.85 ± 0.05 


0.82 ± 0.07 


0.79 ± 0.05 


1.45 ± 0.3 


RMS (m/s) 


5.13 




6.14 


8.00 


6.57 


Stellar "jitter" (m/s) 


2.76 




2.49 


2.02 


5.10 


FAP (from periodogram) 


4.09 X 10"^ 


1.15 X 10"" < 10"" 


8.53 X 10"* 




FAP (from bootstrap) 


< 10-^ 


< 10"* 


< 10"-* 


< 10"* 





Table 6: Derived planet properties. 



^^Parameters include subtraction of linear RV trend in residuals 
■^Evaluated at the time of the first RV point reported for each system 



